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1 Introduction 

The Helmholtz conditions, for the purposes of this paper, are the necessary and sufficient 
conditions for a given system of second-order ordinary differential equations f a (q, q, Q, t) = 
to be of Euler-Lagrange type, that is, for there to exist a Lagrangian A(g, q, t) such 
that 

fa = -rA^)-^- (1) 



dt \dq a J dq a ' 

Here q a are the generalized coordinates (collectively abbreviated to q), q a the corre- 
sponding generalized velocities, and so on. We shall state the conditions shortly. The 
Lagrangian is supposed to be of first order, that is, independent of q and higher-order 
derivative (or more properly jet) coordinates. 

In two recent papers the problem of finding analogous necessary and sufficient conditions 
for a given set of functions f a (q, q, q, t) to take the more general form 



dt \dq a J dq a dq a 

for first-order functions A (a Lagrangian) and D (a dissipation function) has been dis- 
cussed. We shall say that in this case the equations f a = are of Lagrangian form with 
dissipative forces of gradient type. A set of necessary and sufficient conditions is given, 
in terms of standard coordinates, in the fairly recent paper |3j. In a very recent paper [2] 
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a version of the conditions expressed in terms of quasi-velocities, or as the authors call 
them nonholonomic velocities, is obtained. We shall quote the conditions from [3] ex- 
plicitly in Section 2. These conditions are described as generalized Helmholtz conditions 
to distinguish them from the Helmholtz conditions discussed in our opening paragraph, 
which may be called the classical Helmholtz conditions; these must of course comprise a 
special case of the generalized conditions. 

The main purpose of the present paper is to point out that the generalized Helmholtz 
conditions as stated in [3] are not independent: in fact two of them are redundant, in 
that they can be derived from the remaining ones. This we show in Section 2 below. We 
use the same formalism as [3]. Since the version of the generalized conditions obtained 
in [2] is equivalent to that in [3] the same redundancy is present there as well. By taking 
advantage of the improvement in the formulation of the generalized Helmholtz conditions 
that we achieve, we are able to give a shorter and more elegant proof of their sufficiency 
than the one to be found in |3j. 

There are in fact several interesting questions raised by the two papers [21 [3], only one 
of which will be dealt with here. In the third and final section of our paper we give an 
outline of these additional points of interest, which will receive a fuller airing elsewhere. 

We employ the Einstein summation convention throughout. 

We end this introduction with a brief summary of the results about the classical Helmholtz 
conditions that we shall need. 

The classical Helmholtz conditions are that the f a should satisfy 

dfa = dfb_ , . 

dja dh =2 ±(dh\ f) 

dq b dq a dt\dq a J V ' 

9 fa 9f b _ x d { df a df b \ 

dq b dq a ~ 2 dt \dq b dq a ) ' U 
It is a consequence of these (and not an extra condition, as stated in [3]) that 

dq b dq c 

(this follows from the vanishing of the coefficient of q c in condition ([!])). Thus we may 
write f a = g a bQ b + h a , the coefficients being of first order, with g a b = g ba as a result of 
condition ([3]). The Helmholtz conditions can be re-expressed in terms of g ab (assumed 
to be symmetric) and h a , when they reduce to the following three conditions: 

dg a b dg a 



dq c dq 



b 



(6) 



dh a dhb d 

W + W = dTt i9ab) (7) 

dh a _ dhb\ _ d_ ( dha _ dhb\ 

dq b dq a J ~ dt\dq b dq a J ' U 
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where 

That is to say, the f a are the Euler-Lagrange expressions of some first-order Lagrangian 
if and only if f a = g a bq b + h a for some first-order functions g a b and h a such that g a b = 9ba 
and ([S])- ([5]) hold. This reformulation can be found in the book by Santilli [6], for example. 



2 The generalized Helmholtz conditions 

We next turn to the analysis of the generalized Helmholtz conditions. Following the 
notation of [3] we set 

_ dfg df b 1 d ( df b df a 
Tah ~ dq b 8q a 2 dt\dq a dq b 

Sab 2 \dq b + dq a J dt\8q a 

The generalized Helmholtz conditions as given in [3] are that r a b and s a b are of first 
order, and that in addition 

dfa = dfb () 

Qqb Qqa y > 

9s a b 9s ac (-\c\\ 

dr ab _ ds ac _ ds bc , . 

dq c ~ dq b dq a 1 ' 

° = ~dq^ + ~dq^ + ~dq T ' ^ 

Our main concern will be with analyzing conditions (|10|) — f| 1 2 1 ) . which correspond to (2.3e), 
(2.3f) and (2.3g) of [3j; we shall show that conditions (fTUj) and (fT2|) are redundant, being 
consequences of the remaining conditions. 

Our first aim is to understand exactly what it means for r a b and s a b to be of first order, 
bearing in mind condition ([9]) above. From the vanishing of the coefficient of q c in s a b 
we have 

02 f - =0. 

As before f a = g a bQ b + h a , the coefficients being of first order and g a b symmetric. The 
coefficient of q c in r a b is 

1 d ( dg b d .. d , dhb _ dg ad ,. d _ dK 

2 dq c V dq a q dq a dq b Q dq b 

whence 

dgbc = dg ac 

Qqa Qqb 
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The coefficient of q c in s ab , namely 

dg ac dgtc _ 2 dg ab 

Qqb Qga QgC ' 

vanishes as a consequence. The coefficient of q c in r ab is 

dg ac _ dgbc _ i ( d 2 h a d 2 h b \ 
dq b g q a 2\dq b dq c dq a dq c J ' 

an expression which for later convenience we write as p abc ] we must of course have 
Pabc = 0. The remaining terms in r ab and s a b are all of first order, and we have 

dh a _ dhb _ id_ ( dha _ dh b \ 
dq b dq a 2 dt \ dq b dq a ) 

1 ( dha dh\ _ d_, v 

2 { d qb + dq a) dt ySabY, 

compare with ([7|) and ([8]), and also with equations (2.16b) and (2.17c) of [3]. 

The redundancy of condition (fTQj) is a consequence of the vanishing of p abc , as we now 
show. We have 

ds ac _ i / d 2 h a d 2 h c \ _ d_( dg ac \ _ dg ac 
dq b ~ 2 \dq b dq c dq a dq b J dt \ dq b J dq b ' 

using the commutation relation 

" d d_ 

It follows that ds ac /dq b — ds bc /dq a = —p a bc = 0. That is to say, condition (flOl) holds 
automatically as a consequence of the first-order property. Furthermore, p a bc = is 
equivalent to equation (2.17b) of PJ; in other words, the redundancy of (fTUj) is actually 
implicit in [3], though not apparently recognized there. 

Before proceeding to consider condition (|12p we turn aside to make some remarks about 
the classical Helmholtz conditions. The calculations just carried out are essentially the 
same as those which lead to the version of the classical Helmholtz conditions given in 
equations ©-([8]) at the end of the introduction. It is easy to see that in that case 
r ab = s ab = are necessary conditions. This observation, together with the part of the 
argument concerning the vanishing of the coefficients of q and q, leads to the following 



r a b = 
Sab = 



d 



(13) 
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conditions: 



ogac agbc = agah , . 

dq b + Q q a g q c y > 

dh a dh b d 

W + W = 2 dt i9ab) (15) 

(16) 



dg a b 


dg ac 


= 






dq c 


dq b 






d 2 h a 


d 2 h b 




dg ac 


dgbc 


dq b dq c 


dq a dq c 




.dq b 


dq a 


J dK 


dh b \ 


_ d 


(dh a 


dh b 


\dq b 


dq a ) 


~ dt 


\dq b 


dq a 



(17) 
(18) 

These are the conditions quoted in Remark 3 of Section 1 of [3j. However, it is now 
evident that two of them are redundant. Clearly condition (|14p (which is the vanishing 
of the coefficient of q c in s ab ) follows from condition f 1 1 6 [) and the symmetry of g ab - 
Condition (fTT|) is the condition p abc = 0. The second part of the argument above, 
that leading to the relation p abc = ds bc /dq a — ds ac /dq b , shows that in the classical case 
condition (|17p follows from the other conditions. When these two redundant conditions 
are removed we obtain the classical Helmholtz conditions in the form given at the end 
of the introduction. 

These results in the classical case are actually very well known, though not apparently 
to the authors of [3], and have been known for a long time: they are to be found, for 
example, in Santilli's book of 1978 [6] (which is in fact referred to in [3]). For the sake 
of clarity we should point out a difference between the two cases: in the classical case 
condition (fTT|) is completely redundant; in the generalized case it is not redundant, but 
occurs twice in the formulation of the conditions in [3] , once in the requirement that r ab 
should be of first order and once as the condition ds ab /dq c = ds bc /dq a . 

We now return to the generalized conditions, and prove that condition (I12p follows from 
condition ([lip . It will be convenient to write condition (|12p as 

a,b,c 

where ^ a b c stands for the cyclic sum over a, b and c, here and below. As a preliminary 
remark, note that if k abc is symmetric in b and c (say) then ^ a b c k abc = b c ^hac- Now 

dr ab _ d 2 h a d 2 h b 1 d ( d 2 h a d 2 h b 



dq c dq b dq c dq a dq c 2 dt \dq c dq b dq c dq a 
and so by the preliminary remark 

dr ab i d [ v-^ / d 2 h a d 2 h b \\ 1 d ( ^ / d 2 h a d 2 h a 



E 




dq c 2 dt \^ V dq c dq b dq c dq a J j 2 dt 1 V V dq c dq b dq b dq c 

a,b,c ' " 
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On the other hand, using the commutation relation (|13p and the fact that d/dq a and 
d/dt commute it is easy to see that condition (jlip leads to 

1 \ - ( d 2 h a d 2 h a \ _ d 

2 2^ I dq c dq b dq bQ q c )- dt W<*c) ~ U - 
a,b,c ' 



We therefore reach the following proposition, which is stronger than the corresponding 
result in [3]. 

Proposition. The necessary and sufficient conditions for the equations f a (Q,Q,Q,t) = 
to be of Lagrangian form with dissipative forces of gradient type as in |Ij) are that the 
functions r ab and s a b are of first order, that 

dfa = dh 
dq b dq a ' 

and that 

dr ab _ ds ac ds bc 

dq c ~ dq b dq a ' [ ' 

Just as in the classical case we can give an equivalent formulation of these conditions in 
terms of g ab and h a . Bearing in mind that r ab and s ab being of first order are essential 
hypotheses, we find that the following conditions are equivalent to those given in the 
proposition above: f a = g a bQ b + ha with g ab symmetric, where g ab , h a are of first order 
and further satisfy 



dg a b dg a 



dq c dq b 

d 2 K d 2 h b 







I be 



E 

a,b,c 



dq b dq c dq a dq c \ dq b dq a 
d 2 h a d 2 h a 



dq b dq c dq c dq l 



0. 



The first of these is one of the classical conditions. The second is the condition p abc = 0, 
which holds in the classical case as we have shown. The third is just condition f)19f) 
above expressed in terms of g ab and h a (or as it turns out, in terms of h a alone), and 
r a b = s ab = in the classical case. It is evident therefore that the conditions above are 
indeed a generalization of those for the classical case. 

We end this section by giving an alternative proof of the sufficiency of the generalized 
Helmholtz conditions, based on this formulation of them, which is shorter and in our 
view more elegant than the proof in [3] (necessity is an easy if tedious calculation) . 

We note first that if g ab is symmetric and satisfies dg ab /dq c = dg ac /dq b then 

_ d 2 K 

9ab ~ 8q a dq b 
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for some function K = K(q,q,t) (a well-known result, which also appears in [3]). Of 
course K is not determined by this relation; in fact if A = K + P a q a + Q, where P a and 
Q are any functions of q and t, then A has the same Hessian as K (the same g a b, in other 
words). Our aim is to choose P a and Q so that the given equations are of Lagrangian 
form with dissipative forces of gradient type as in ([2]), with Lagrangian A, assuming that 
the generalized Helmholtz conditions above hold. In fact we won't need to consider Q 
because it can be absorbed: if A is a Lagrangian and D a dissipation function for some 
functions f a , so are A + Q and D + q a dQ/dq a . We shall therefore take Q = below. 

Let E a be the Euler-Lagrange expressions of K. Then E a = g a bq b + k a for some first- 
order k a , by construction, so f a — E a = h a — k a = K a say, where n a is also of first order. 
Moreover, f a satisfies the generalized Helmholtz conditions by assumption, and E a does 
so by construction (it satisfies the classical conditions after all), whence K a satisfies 



d 2 K a d 2 K b 



(20) 



QqbQqC dq a dq c 
V ( _ 9 2 Ka \ = 

f^ c \dq b dq c dq c dq h ) ' 1 ' 

Let us set dK a /dq b — dn^/dcf 1 = R a t,. Then by (f20|) R a i, is independent of q, and by (f2Tj) 

E dRgb _ q 
dq c ~ 

a,b,c 

There are therefore functions P a (q,t) such that 

dP a dP b \ dn a dK b 



Rab — 2 

which is to say that if we set 



Qqb Qqa J Qqb Qqd 



dKg _ fOPa _ dPh 

nab dq b \dq b 8q a 

then ir ha = T^ab- Moreover, 

dTT ab d 2 K a dTT ac 



dq c dq b dq c dq b 
It follows (just as is the case for g a b) that there is a first-order function D' such that 

8 2 D' 

nob ~ WW' 

from which we obtain 

_(dPa_dPb\ b + W +S 
a ~\dq» dq-) q + dq« + ba 
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where S a is independent of q. Now take A = K + P a q a - Denoting the Euler-Lagrange 
expressions of K by E a as before, the Euler-Lagrange expressions for A are 



dP a dP b \ , 6 , dP a „ , 3D' „ , <9P a 



W~WJ 4 + ~df = Ea + Ka ~W~ Sa + ot 



Thus, putting 



d fdA\ dA dD 



we get 

f - = j 

dt \dq a J dq a dq a 

as required. 

This method of proof works equally well in the classical case. The proof is constructive, 
in the same sense that the one in [S] is, in either case. It is particularly well adapted 
to the familiar situation in which g a b is independent of q, when one can take the kinetic 
energy \g a bq a q b for K. 



3 Concluding remarks 

We wish to make four remarks in conclusion. 

The first remark concerns the nature of conditions (|10p - (]12p on the derivatives of r a b 
and s a b, as originally expressed in [3] (that is, ignoring the question of dependence). In 
particular, bearing in mind the fact that r ao is skew in its indices, the condition 

Q q c 

a,b,c 

is suggestive: if perchance the were functions of the q alone this would have a nat- 
ural interpretation in terms of the exterior calculus, being the condition for the 2-form 
r abdq a A dq b to be closed, that is, to satisfy d(r a bdq a A dq b ) = 0. This point is made, 
in somewhat different terms, in [3] (and we appealed to the same general result in our 
proof of sufficiency of the generalized Helmholtz conditions in Section 2). The authors of 
[3J go on to say, however, that the condition above 'can be interpreted as the vanishing 
curvature of a symplectic space', which seems to us not to be entirely convincing. In 
fact it is possible to interpret the three conditions (|10 p - (|12p collectively as signifying the 
vanishing of a certain exterior derivative of a certain 2-form on the space of coordinates 
9> Qy 9) • • • j a 2-form whose coefficients involve both r a b and s Q &. This interpretation re- 
ally arises from seeing the problem in the context of the so-called variational bicomplex 
(see [7] for a recent review). 

Secondly, we contend that the problem we are dealing with should really be regarded 
as one about (second-order) dynamical systems. The point is that a dynamical system 
may be represented as a system of differential equations in many different coordinate 
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formulations; the question of real interest is whether there is some representation of it 
which takes the form of an Euler-Lagrange system with dissipation, not just whether a 
given representation of it takes that form. Of course this point applies equally, mutatis 
mutandis, to the case in which there is no dissipation. Now the Helmholtz conditions 
as discussed in [2JE], in both the classical and the generalized versions, suffer from the 
disadvantage that they are conditions for a given system of differential equations to be of 
Euler-Lagrange type. There is, however, an alternative approach to the problem which 
does deal with dynamical systems rather than equations, at least in the case in which 
the system can be expressed in normal form q a = F a (q,q,t). In this approach one asks 
(in the absence of dissipation) for conditions for the existence of a so-called multiplier, a 
non-singular matrix with elements g a b, such that g a b{Q b — F b ) takes the Euler-Lagrange 
form (so that in particular when the conditions are satisfied g a b will be the Hessian of the 
Lag rangian with respect to the velocity variables). The basic idea is to put h a — — gab-F 
in the conditions at the end of the introduction, and regard the results as a system of 
partial differential equations equations for g a ^ with F a known. The seminal paper in 
this approach is Douglas's of 1941 [I], which analyses in great detail the case of two 
degrees of freedom. For a recent review of developments since then see Sections 5 and 6 
of [4] and references therein. One can in fact also formulate conditions on a multiplier 
for a second-order dynamical system, expressible in normal form, to be representable as 
equations of Lagrangian form with dissipative forces of gradient type; these generalize the 
known results for representation in Lagrangian form without dissipation in an interesting 
way. 

The new ingredient in [2J, by comparison with [3], is the expression of the generalized 
Helmholtz conditions in terms of quasi- velocities. As presented in the paper this is quite 
a long-drawn-out procedure, because in effect the conditions are rederived from scratch. 
Our third remark is that in principle this should be unnecessary: a truly satisfactory 
formulation of the conditions should be tensorial, in the sense of being independent of 
a choice of coordinates (and of course quasi-velocities are just a certain type of velocity 
coordinates). The approach described in the previous paragraph leads to conditions 
which have this desirable property. 

Fourthly and finally, there is the question of whether generalized Helmholtz conditions 
can be derived for other kinds of "generalized force" terms than dD/dq a . One important 
case is that in which such a term is of gyroscopic type. We have obtained such conditions 
in this case, again using the approach discussed in our second remark above. 

These points are discussed in full detail in a recently written paper [5] . 
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